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Josephson junctions with a ferromagnetic metal weak link reveal a very strong decrease of the critical current 
compared to a normal metal weak link. We demonstrate that in the ballistic regime the presence of a small 
region with a non-collinear magnetization near the center of a ferromagnetic weak link restores the critical 
current inherent to the normal metal. The above effect can be stimulated by additional electrical bias of the 
magnetic gate which induces a local electron depletion of ferromagnetic barrier. The underlying physics of 
the effect is the interference phenomena due to the magnetic scattering of the Cooper pair, which reverses its 
total momentum in the ferromagnet and thus compensates the phase gain before and after the spin-reversed 
scattering. In contrast with the widely discussed triplet long ranged proximity effect we elucidate a new 
singlet long ranged proximity effect. This phenomenon opens a way to easily control the properties of SFS 
junctions and inversely to manipulate the magnetic moment via the Josephson current. 

Mesoscopic properties of nanometer sized conductors are strongly affected by the injection of correlated 
electrons (Cooper pairs) from superconducting electrodes. While propagating through normal metal such 
pairs of electrons are able to preserve their superconducting correlation on mesoscopic lengths providing 
the superconducting current flow through SNS (superconductor-normal metal-superconductor) weak links 1,2 . 

Time reversal of electronic states forming a Cooper pair is an important component of the above correlation. 
Absence of superconducting pairing interaction in a N part of the superconducting SNS device opens a possibility 
of easy external manipulation of the spin structure of the propagating Cooper pairs. This offers the means for 
spintronic manipulation of superconducting weak links. The most efficient "tailoring" of Cooper pairs can be 
achieved by external in-homogeneity located on a submicron length scale, which is set by superconducting 
coherence length, determining the scale of spatial correlation of paired electrons. Here we demonstrate that a 
tip of the magnetic exchange force microscope (MExFM) 3-4 , which induces localized in space magnetic exchange 
fields can play the role of such a local probe for a spin state of superconducting Cooper pairs. Existing experi- 
mental evidences of the externally induced exchange fields h in metals of the order of few 5,6 or even few tens 7 
millielectronvolts place the electronic coupling to such field in a range of energies comparable with (or even 
exceeding) superconducting energy scale ~ 1 me V. We will show that the effect of exchange induced gating of 
Cooper pairs leads to a new phenomenon - stimulation of long- range singlet superconductivity in SFS (super- 
conductor-ferromagnet-superconductor) weak links. 

Magnetic exchange interaction in ferromagnetic metals lifts a degeneracy with respect to spin orientation of the 
electrons, forming a Cooper pair. This leads to different de-Broglie wavelengths of electrons at Fermi surface for 
spin-up and spin-down orientation and produces a modulation of the Cooper pair wavefunction while prop- 
agating along the ferromagnet 8 . As a result an oscillatory damping of the superconducting ordering is known to 
appear when a ferromagnetic ordering occurs in a normal metal link connecting two S electrodes 9 . This phe- 
nomenon provides the basis of the n -junction realization 1012 . Considering the quantum mechanics of quasipar- 
ticle excitations the exchange field leads to phase difference y ~ L/^f, gained between the electron- and hole- like 
parts of the total wave function along a path of the length L 13,14 , where Cji = KVr/2h is a characteristic length 
determined by the exchange field (V F is the Fermi velocity). 
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Figure 1 | (A) A schematic picture of singlet Cooper pair scattering with no 
spin-flop transition of electrons. The spin arrangement of a pair/+ (q = kj 
— k|) does not change with respect to the exchange field: q' = q. (B) A 
schematic picture of "exchange" singlet Cooper pair scattering with spin- 
flop transition of electrons: /+—»/_. The spin arrangement of a pair /_ was 
changed with respect to the exchange field: q' = k t k T = q. 
The scattering domain d 2 is shown by grey color. The symbols /+ indicate 
the Cooper pairs with zero spin projection and a reversed spin 
arrangement (see Methods). 

Measurable quantities should be calculated as superpositions of 
fast oscillating contributions e' y from different trajectories and, thus, 
rapidly vanish with the increasing distance from the SF boundary. It 
should be noted though, that a simple domain structure consisting of 
two F layers with opposite orientations of exchange field cancels the 
phase gain y and suppresses the destructive effect of an exchange 
field 1315 . It is clear therefore that by affecting the spin structure of 
Cooper pairs, one influences the strength and spatial distribution of 
the proximity effect induced by Cooper pair penetration inside a 
ferromagnetic metal. It was suggested before that the Cooper pairs 
of electrons with aligned spins (with equal-spin pairs) can be formed 
by spatial non-homogeneous magnetization 1617 . Since they bind elec- 
trons with exactly the same de-Broglie wave length, these triplet 
Cooper pairs do not dephase, thereby leading to long-range prox- 
imity effect. Singlet Cooper pairs still become "filtered out" from 
spatial transfer of superconducting phase coherence due to a strong 
dephasing effect, occurring in long (as compared with magnetic 
coherence length ^) SFS weak links. To observe such a long ranged 
triplet superconducting current, the SFS junctions with a composite 
F layers comprising three non-collinear domains, were suggested 
theoretically 18-20 and realized in recent experiments 21,22 . In such a 
case the triplet component is generated by a thin ferromagnetic 
domain, located between superconducting lead and a thick central 
non-collinear domain. The long ranged Josephson current results 
from the interference between these triplet components, generated 
by the left and right superconducting leads. 

Here we suggest a new way of manipulating the Cooper pairs flow 
through a ferromagnet, which consists in using a well controlled tip/ 
probe along the supercurrent flow. This corresponds to the case of 
the composite F layer with a thin domain, located near the center of 
the junction. In contrast with the situation analyzed in Refs. 18-20, 
the triplet superconducting current is absent for this setup but 



becomes possible a long-ranged singlet proximity effect. As we will 
show the field generated by the probe induces a special scattering of 
Cooper pairs which corresponds to exchange spins of two electrons 
forming a pair. A schematic picture of such scattering is shown in 
Fig. 1. As we have previously mentioned two electrons forming a 
singlet Cooper pair have a non-zero total momentum fiq = fikf — fik| 
due to the ferromagnetic exchange splitting of the spin subbands 
(The modulus of the Fermi wave-vector for electrons with a spin 
polarized along the field is larger |kj| > |kjj and |q| ~ The 
electrons in a singlet Cooper pair reach the scattering center (spin 
exchanger) and scatter their spin so that the new total momentum of 
the Cooper pair hq' is either unchanged Hq' = fiq (Fig. 1 A) or 
reversed Sq'= — hq (Fig. IB) (see the discussion in Methods). In 
the first case the spin arrangement of a singlet Cooper pair has not 
changed with respect to the exchange field and there remains a total 
phase gain y ~ {d\ + d 3 )/^i, which results in a strong suppression of 
proximity due to the destructive trajectory interference. In the sec- 
ond case the scattered Cooper pair has a reversed spin arrangement 
and the total phase gain is y ~ (d 1 — d 3 )/^i,. As a result, at a sym- 
metric position of the scatterer (di = d 3 ) the total phase gain for a 
singlet Cooper pair should be cancelled (y = 0) and the long range 
singlet superconducting proximity in SFS link becomes possible. 

To be more precise, we consider the Josephson transport through a 
normal ballistic nanowire (NW) in contact with a ferromagnetic 
insulator (FI). The FI turns the NW into an effective ferromagnet 
with an exchange field h. The schematic picture of the SFS device is 
presented in Fig. 2A. The total length of the constriction d = d\ + d 2 
+ d 3 is assumed to be large compared to the magnetic coherence 
length C), =hVp/2h : d2>C(,- For simplicity we restrict ourselves to 
the case of a short junction with <i<Cf„, where = HVp/T c is the 
coherence length of normal metal (T c is the critical temperature of 
the S layer). The magnetic tip is assumed to bring on localized in 
space magnetic exchange field inhomogeneity which we model by a 
stepwise profile: 

( hz a , in domains d\,dx 

Hz) = \ , \ (1) 

[ h (z 0 cos a + x 0 sin a) , in domain d 2 , 

where a is the angle of the exchange field rotation in the central 
domain d 2 (see Fig. 2B). 

Results 

The current-phase relation of SFS Josephson junction is determined 
by the quasiclassical relation 14,15 

((n, n f )cos ny) 
' <(n,n f )) ' 



sin n<p - 



(2) 



where n is the unit vector normal to the junction plane, n F is the unit 
vector along the trajectory, and a„ are the coefficients of the Fourier 
expansion for the current-phase relation for superconductor-nor- 
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Figure 2 | (A) The schematic sketch of the SFS constriction under consideration: normal metal nanowire (NW) in contact with a superconductor (S) 
and a ferromagnetic insulator (FI). (B) Equivalent SFS Josephson junction containing three ferromagnetic layers (domains) with a stepwise 
profile of the exchange field (1). Linear quasiparticle trajectory is shown by the red dashed line. 
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mal metal junction of the same geometry. The angular brackets 
denote the averaging over different quasiclassical trajectories char- 
acterized by a given angle 0 and a certain starting point at the super- 
conductor surface, and for 3D constriction looks as 



((n, rif)cos(ny)) 
<(n,n P )} 



It/2 

= 2 1 dO sin 0 cos 0 cos(ny) 



(3) 



where cos 0 = (n, n F ). At temperatures T close to T c the current- 
phase relation (2) is sinusoidal, and the coefficient a t is determined 
by the following simple relations 15 : 



«i ■ 



8h 



-N 



(4) 



Here A is the temperature dependent superconducting gap, The 
factor N is determined by the number of transverse modes in the 

junction: N = s 0 ~ 1 j ds j drip(rip, n) ~S/s 0 , where S is the junction 

cross-section area, and s^ 1 = (kp /2tz) 2 , where k F is the Fermi 
momentum. 

The phase gain y can be conveniently determined from the 
Eilenberger — type equations 23 if we use a standard parametrization 24 
of the anomalous quasiclassical Green function/ =f s + tt o, where f s 
(f ( ) singlet (triplet) parts of the function, respectively, and a is a Pauli 
matrix vector in the spin space. The functions/, f, satisfy the linear- 
ized Eilenberger equations 18 written for zero Matsubara frequencies 

- ihV F d s f s + 2hf , =0, - ihV F d s f, + 2f s h = 0, (5) 

with the boundary conditions/(s = s L ) = 1, f f (s = s L ) = 0 at the left 
superconducting electrode (for simplicity we consider the case of the 
absence of the barriers at the interfaces). The phase gain y along the 
trajectory in the equivalent SFS junction (Fig. 2B) is determined by 
the singlet part of the anomalous quasiclassical Green function/(s = 
s R ) = cos y taken at the right superconducting electrode 15 . Solving the 
equations (5) for the stepwise profile of the exchange field (1) we find 
(see Methods): 



cos)» = cos £> 2 cos(<5i + <5 3 ) — cos a. sin 8 2 sin(<5i +<53) 
— sin 2 asin(5i sin 63(1 —cos 82), 



(6) 



where cos 0 = (n, n f ) and <5, = dj^ h cos 0 (i = 1, 2, 3). Averaging the 
expression (6) over the trajectory direction 0 and neglecting the 
terms proportional to {^/d<Cl, which decrease just as for the case 
of homogeneous ballistic 3D SFS junction, one arrives at the follow- 
ing long-range (LR) contribution: 



1 \LR 
(cosy) =• 



1 



sin a(l — cos <5 2 )cos 28, 



(7) 



where 8 Z = z 0 /^h cos 0 and z 0 = {di — d 3 )/2 is the shift of the central 
domain with respect to the weak link center. So, the long-range 
component of the Josephson current at the first harmonic is deter- 
mined by the relation: 



T LR 



sin <p, 



n/2 

T™=2j dO sin0 cos 0{cos) 



(8) 



For a very thin central domain c? 2 <C£j, in the center of the NW (z 0 = 
0) one can easily estimate from (7, 8) the critical current of the SFS 
junction 



LR\, 



J 0 . 2 

— sin a 



In 



(9) 



where J 0 = (eT c N/8h) {MT c f - is the critical current of the SNS 
junction for zero exchange field (y = 0). We see that the long-ranged 



critical current reaches the maximum at a = n/2 and grows with the 
increase of d 2 up to d 2 ~ C/,. Our numerical calculations show that it 
is maximum for c2 2 — 2.5 and may reach ~Q.7I 0 . Certainly, the 
above long-range effect in the first harmonic describes the properties 
of the SFS constriction if contribution of higher harmonics in the 
current-phase relation (2) is negligible. We present the analysis of 
the second harmonic effect in the Supplementary information. 



Discussion 

Figure 3 shows the dependence of the maximal Josephson current 
I^ R = aiT^ R on the thickness d 2 of the 90° domain (a = n/2) for 
different values of the shift of the domain z 0 . Naturally, when the 
thickness of the central domain d 2 goes to zero, the long-range effect 
disappears. We can see from Fig. 3, that the long-range component 
of the Josephson current I LR coincides approximately with the total 
supercurrent across the junction (2) until the outer domains are long 
enough: d\, d^^^. The amplitude of I LR depends nonmonotonically 
on the size of the central domain d 2 and has the first maximum at 
d 2 — 2.5^1,. Interestingly, that the long-range contribution generates 
a re— junction (I LR is negative for zero shift of the domain z 0 = 0). 
With the shift of the central domain the junction can be switched 
from 7t— to 0— state. Figure 4 shows the dependences of the maximal 
Josephson current I c on the position of the central domain z 0 for 
different values of the rotation angle a. We may see that the critical 
current is very sensitive to the position of the central domain and the 
first zero of 1\ occurs already at z 0 — 0.5^. 

Here we considered a simple step-like model of magnetization 
distribution in F layer. For the very thin central domain d 2 <^i^, 
we may easily estimate the long ranged contribution (cos y) LR ~ 
(d 2 /£h) 2 (h x /h) 2 which is in accordance with the expression (7). For 
a general profile of the magnetization it may be convenient to use the 
transfer-matrix method (see the Supplementary information). The 
smooth (on the scale ^) profile of the magnetization decreases the 
long ranged effect and the proposed mechanism occurs to be most 
efficient for d 2 ~ C^. 

Note that the considered phenomenon should generate the oscil- 
lating potential profile for the magnetic tip U(z 0 ) ~ — I c (z 0 ) cos (p 
which depends on superconducting phase difference across the junc- 
tion. This opens an interesting possibility to couple the Josephson 
current oscillations with mechanical modes of the tip. On the other 
hand the same effect can produce a change of the orientation of the 
magnetic moment. Inversely, the precession of the magnetic moment 
shall modulate the critical current of the junction and provides a 
direct coupling between the superconducting current and the mag- 



0.5 




Figure 3 | The dependence of maximal Josephson current = maxf/Jon 
the thickness d 2 of the 90 ° domain (or = jtI2) for different values of the 
shift of the domain Zq: Zq = 0 - red solid line; Zq = g h - blue dashed line; Zq 
= 2ft - green dash-dotted line. Symbols + show the long-range part of the 
supercurrent. Dotted line shows the value of I c in absence of domain d 2 (a 
= 0). We have set T = 0.9T C ; d = 20o, [I 0 = (eT c N/8h) (D/T c ) 2 ]. 
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formalism. For this, we need to solve the linearized Eilenberger equations written for 
zero Matsubara frequencies 



- iW P dJ s + 2hf , - 0, - iHV F d s f, + 2/ s h - 0, 



(10) 



Figure 4 | The dependence of maximal Josephson current = max{/]) on 
the shift of the central domain Zo for different values of the angle a:a=nl 
2 - red solid line; a = nli - blue dashed line. Symbols + show the long- 
range part of the supercurrent. We have set T= 0.97^ d= 20^; d 2 = 2.5c;„ 
U 0 = (eT c N/8h) (A/T c ) 2 ]. 

netic moment in the weak link similar to the situation discussed in 
Refs. 25,26. 

The magnetically tunable long-range SFS proximity effect sug- 
gested above has a potential to be an important feature of carbon- 
based superconducting weak links. Graphene sheets and carbon 
nanotubes are reported to offer a ballistic propagation for electrons 
on a micrometer length scale 27-28 . This fact together with appearing 
reports on a gate tunable magnetism in graphene 29,30 makes all ingre- 
dients of the present theory achievable in experiment. Another pos- 
sibility is to use the indium antimonide (InSb) nanowires as a 
superconducting weak link. The indium antimonide nanowires, 
recently used in the experiments to reveal the signature of 
Majorana fermions 31 , demonstrated a very high ^-factor (j~50). 
Anomalously large g-factor reported in such wires offers the pos- 
sibility to "mimic" a ferromagnetic spin-splitting effect of the order 
of 10 K by simply applying an external magnetic field of the order of 
0.1 Tesla, and then making such nanowire a suitable candidate for a 
weak link to observe the discussed phenomena. Note that in contrast 
to the experiments 31 the magnetic field should be applied along the 
spin-orbit field axis to avoid the interference with the spin-orbit 
effect. 

It should be noted, that new additional functionality of the con- 
sidered device can be achieved by electric biasing of the magnetic 
gate 32,33 . In weakly doped ferromagnetic barriers, such bias ( V g ) alters 
both the charge carrier concentration and the Fermi velocity. 
Choosing a polarity of electric gating one can create a depletion 
region beneath the tip. As a result, both the Fermi velocity V F and 
the exchange length £ h = hVp/2h decrease in the spatial region of the 
domain d 2 , and the key parameter responsible for the magnetic 
exchange scattering 5 2 = d 2 lc,h grows. For thin domains (c^-Cl) 
the critical current I c ~d 2 increases with the gate voltage V g , and 
the local depletion of F barrier should result in the stimulation of 
the superconductivity. This nontrivial interplay between electric and 
magnetic gating effects can be used to control singlet Josephson 
current through ferromagnetic nanowires. 

To summarize, we studied the interference phenomena originated 
by the spin-exchange scattering in ferromagnetic ballistic weak link 
and demonstrated that they provide an efficient way to control the 
Josephson current and to couple it with a magnetic moment. 

Methods 

A. Transfer-matrix formalism for Eilenberger equations. To consider the 
Josephson transport through ferromagnetic layer with a non-collinear 
magnetizations M and exchange field h it is convenient to utilize the transfer-matrix 



for the case when the quantization axis is taken arbitrarily in the ferromagnetic layer 
of a thickness d. We assume that a quasiclassical trajectory s and exchange field h = h 
(zo sin a + Xo cos a) lie in the plane (x, z), as shown in Fig. 5. The trajectory is 
characterized by a given angle 6 with respect to the z-axis. The triplet part f ( of the 
anomalous quasiclassical Green function/ —f s + i t d consists of two nonzero 
components and can be written as f f = / (JC Xo + f tz z 0 . Defining the transfer-matrix 
T^(d,6) that relates the components of the Green function/(s) at the left (5 = 0) and 
right (5 — SJ = d/cos 8) boundaries of the F layer, 

//.(*) \ 

}(s d )= Ms*) =r»(d, 

WW 

we get the following expression: 

/ cos(qs d ) 




(ID 



-z'cosasin(gsj) 



(12) 



f a (d,6) — — i cosa sin(gs</) sin 2 a + cos 2 a CQsiqsj) 
\ — isinasin(^5j) sinacosafcos^sj) — 1) 

— fsina sin(gsd) \ 
sin a cos v.(cos(qsd) — 1) 
cos 2 oc + sin 2 a cos(qsd) j 

where q = l/£ h — 2h/HV F . 

In order to elucidate the peculiarities of the Cooper pairs scattering with a spin-flop 
transition of electrons it is convenient to introduce the new functions f ± = f s ± f tz 
which describes the pairs with zero spin projection and a reversed spin arrangement. 
The transfer-matrix T ri (d,9) can be drastically simplified if the direction of the 
exchange field coincides with a spin quantisation axis z. In this case, a = 0 and 
f± ( s d) — e + " ,Sd f± (£>)■> ftx( s d) ~_f tje (0). Calculating the superconducting current at the 
right electrode S# we readily see that it results from the interference with the singlet 
component coming from the left electrode _^(5^) — (f+isj) + /_{s^))/2 (triplet com- 
ponents are irrelevant because the right electrode provides only the singlet com- 
ponent). The oscillating factors e + "* Sd in f±(s ( j) produce, after the averaging over the 
trajectories directions (angle 9), a strong damping of the critical current compared to 
the normal metal (where these factors are absent). 

Now we may easily understand the mechanism of the long-ranged proximity effect. 
Indeed after coming through the first F layer an extra phase factor appears inf± 
functions (see Fig. 2):f+ (sd 1 ) — e +ll,Sdl f s (0). In the absence the middle layer, the/ ± at 
the right electrode would be/ ± (s di + s di ) - e + " J ( s "i +Sd ^f s (0) and the oscillating 
factors will strongly damp critical current. The additional non-collinear middle layer 
d 2 will mix up the components/)- and/_ - see the matrix (12) and, for example, 
/+ (s^ +Sd 2 ) in addition to e~"* Sii i_/i(0) component will have a e+'^'^O) contri- 
bution, i. e.f + (s,j, J rSd 2 )=a e~ leJSdl + fee + " ?id| . In fact, namely this mechanism is 
schematically presented in the Fig. 1(b). Then the resulting/ + function at the right 

electrode should be/+(s dl +s d2 +s di )=a +s "^ +be +iq ( s * and for d x = 

d 3 the oscillating factor at the second term vanishes. This means the emergence of the 
long-ranged singlet proximity effect discussed in the present report. Note that the 




Figure 5 | Quasiclassical trajectory s through homogeneous 
ferromagnetic layer of a thickness d with an arbitrary direction of the 
exchange field h. 
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additional noncollinear F layer d 2 may strongly increase the critical current, provided 
that it is placed at the center of the structure. 

The transfer-matrix method is very convenient for the calculation of the Josephson 
transport through the SFS junction containing three ferromagnetic layers with a 
stepwise profile of exchange field. For this geometry shown in Fig. 2B of the Letter, the 
anomalous quasiclassical Green function /(s R ) — (f s (sR),f tz (s R ),f lx (s R )) at the right 
superconducting electrode (s = s R — d/cos 6) can be easily expressed via the boundary 
conditions /(O) — (1,0,0) at the left superconducting electrode (5 = 0) as follows: 



f( SR ) = r o (4,0)r a (4,0)To(*.0)/(°). 



(13) 



where d = d x + d 2 + d 3 is the total thickness of the ferromagnetic barrier. As a result, 
the singlet part^(s^) responsible for the Josephson current through the junction, can 
be written in the form (6). 
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